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Preface 

The  purpose  of  this  study  was  to  obtain  and  analyze  the 
electromagnetic  shadowing  function  for  a  random  rough 
surface  described  by  a  true  exponential  probability  density 
function.  It  has  been  determined  that  large  errors  may 
occur  if  shadowing  is  neglected  when  computing  scattering 
coefficients.  This  report  developed  and  compared 
electromagnetic  shadowing  functions  for  three  different 
probability  density  functions. 

This  project  employed  a  unique  blend  of 
electromagnetic  scattering,  probability  and  advanced 
calculus  principles. 

I  would  like  to  thank  my  advisor,  Dr.  Vital  Pyati,  who 
derived  the  true  exponential  probabi lity  density  function 
used  in  this  thesis,  for  his  timely  assistance  in  the 
development  of  this  thesis.  A  very  special  thanks  goes  to 
Lt  Col  William  Baker.  Without  his  guidance  and  help  with 
the  math,  this  thesis  could  not  have  been  accomplished. 

Also,  I  would  like  to  thank  Dr.  Robert  J.  Papa  of  Rome  Air 
Development  Center  ( RADC ) ,  a  noted  author  in  the  field  of 
rough  surface  scattering,  for  his  valued  inputs;  John,  Pam 
and  the  rest  of  the  outstanding  staff  of  the  AF1T  Library 
for  their  never  ending  support;  and  Capt  Thomas  Roseman  for 
his  can  do  approach  to  problems  posed  in  our  discussions  of 
our  theses. 
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Abstract 

The  purpose  of  this  study  was  to  develop  a  shadowing 
function  for  a  rough  surface  model  that  was  based  upon  an 
exponential  probability  density  function  (pdf).  This  study 
had  four  critical  parts: 

(1)  Derive  the  characteristic  function  and  the  joint  slope 
probability  density  function  for  a  second-order  exponential 
pdf.  (2)  Write  a  computer  program  to  numerically  compute 
the  shadowing  function  for  the  second-order  exponential  pdf 
surface  model.  (3)  Compare  the  shadowing  function  obtained 
with  an  'exponential  like’  and  Gaussian  pdf.  (4)  Also, 
investigate  a  discrepancy  in  the  mean-square  total  slope 


The  study  found  that  the  Shadowing  function  for  the 
exponential  pdf  had  a  lower  value  than  the  shadowing 
function  for  the  'exponential  like’  and  gaussian  pdfs.  The 
shadowing  functions  for  the  'exponential  like’  and  Gaussian 
pdf  were  found  to  have  about  the  same  value.  This  study 
also  found  that  the  shadowing  function  is  sensitive  to  the 
value  of  the  mean-square  slope  term  used. 


ELECTROMAGNETIC  SHADOWING  BY  AN 


EXPONENTIALLY  DISTRIBUTED 
ROUGH  SURFACE 

I ■  I ntroduc t i on 

In  recent  years,  the  performance  of  a  complex  and 
versatile  radar  system  has  become  more  sensitive  to  its 
operational  environment.  This  has  created  a  need  for 
accurate  mathematical  models  of  the  environment  in  which 
radar  systems  are  used.  For  example,  a  model  for  electro¬ 
magnetic  scattering  by  a  randomly  rough  surface  should 
account  for  multipath  effects  and  shadowing  of  the  incident 
as  well  as  the  reflected  wave.  The  rough  surface  models 
must,  of  course,  be  mathematically  sound  and  preferably  in 
agreement  with  measured  data. 

Early  models  used  a  Gaussian  random  probability  density 
function  (pdf)  because  of  their  analytical  convenience 
[9:1274].  However,  different  surface  models  are  needed  for 
different  magnitudes  of  surf  ace  roughness.  For  example, 
Gaussian  random  surface  models  do  not  adequately  model  sea- 
ice  fi.elds  and  some  kinds  of  rough  terrain  because  of  their 
sharp  peaks  and  steep  slopes.  To  accurately  model  these 
surfaces,  a  joint  non-Gauss  ion  random  probability  density 
function  is  preferred  [10:788]. 
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In  this  connection,  Barrick  introduced  the  following 
exponential  joint  surface-height  probability  density 
function  to  model  these  surfaces. 


Here,  K  and  <  ,  which  are  functions  of  x  and  y,  denote  two 

heights  of  the  rough  surface  at  points  1  and  2  which  are 
separated  by  a  distance  r.  P(r)  is  the  surface-height 
correlation  function  and  h2  is  the  mean-square  roughness 
height.  Also,  it  is  assumed  that  P(r)  does  not  contain  a 
steady  or  periodic  component  [16:720,721]. 

Pyati  stated  that  Eq  (1)  does  not  meet  one  of  the 
statistical  requirements  of  a  pdf  and  is  therefore,  math¬ 
ematically  incorrect.  Hence,  Pyati  provided  the  following 
second-order  exponential  pdf  that  does  meet  all  the 
statistical  conditions  [14:8]. 
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Here,  sgn!C)  is  the  signum  function  which  is  equal  to  plus 
one  for  positive  C  and  minus  one  if  <  is  negative.  Ifc  is  a 
modified  Bessel  function  of  the  first  kind  of  order  k.  P  is 
P(r),  which  is  the  surface-height  correlation  function  and 
h2  is  the  mean-square  roughness  height. 


Prob 1 em 

The  purpose  of  this  study  is  to  investigate  shadowing 
using  Pyati's  second-order  exponential  pdf.  Here,  shadowing 
is  defined  as,  "the  screening  of  parts  of  the  surface  by 
other  parts  interrupting  the  illuminating  radiation" 

[4:384).  The  results  of  this  investigation  will  be  compared 
with  those  obtained  by  Brown  when  he  used  Barrick’s  'expo¬ 
nential  like’  pdf,  Fq  (1)  [10:788-790].  It  may  be  noted 


that  among  others,  Brown  has  used  Barrick ’ s  Eq  (1)  in  his 
investigation . 


Assumptions 

It  will  be  assumed  that  the  incident  wave  striking  the 
surface  is  plane.  This  assumption  is  satisfied  by  placing 
the  receiving  and  transmitting  antennas  sufficiently  far 
from  the  rough  surface.  The  wavelength  of  the  incident  wave 
will  be  assumed  to  be  much  shorter  than  the  dimensions  of 
the  illuminated  surface  area.  Also,  the  randomly  rough  sur¬ 
face  will  be  considered  to  be  a  perfect  conductor. 

Scope 

Although  the  problem  of  electromagnetic  scattering  and 
shadowing  is  very  extensive,  this  i n ves t i ga t i on  will  only 
consider  the  study  of  Pyati’s  second-order  exponential  pdf 
following  Brown’s  method.  A  comparison  between  the  shadow¬ 
ing  results  when  using  Barrick 1 s  and  Pyati’s  second-order 
pdfs  will  be  made. 

Shadow  j  ng  Overv i ew 

Scattering  coefficients  of  rough  surfaces  are  prone  to 
large  errors  if  shadowing  is  neglected  [4:384].  To  correct 
for  these  errors,  Beckmann  used  a  shadowing  function  that 
was  equal  to  one  for  illuminated  areas  and  equal  to  zero  for 
shadowed  areas  [  1:384]  .  However,  because  of  its  neglect  of 
the  spectrum  of  shadows  caused  by  diffraction,  Beckmann’s 
shadowing  function  is  only  valid  in  the  high  frequency  limit 
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where  the  wavelength  approaches  zero,  which  is  the  same 
restriction  geometrical  optics  has  [4:388].  Beckmann 
concluded  that  fields  scattered  by  a  rough  surface  are 
strongly  affected  by  shadowing,  especially  for  large  angles 
of  incidence . 

Expanding  on  Beckmann’s  work,  Smith  developed  a  shadow¬ 
ing  function  based  upon  the  geometry  of  self  shadowing. 
Smith’s  results  agreed  closely  with  a  simulation  of  a  Gaus¬ 
sian  random  surface  by  Brockelman  and  Hagfors  [19:668]. 

Brown  simplified  Smith’s  shadowing  function  by  assuming  that 
the  height  and  slopes  of  the  surface  were  independent  random 
variables.  Using  this  assumption,  Brown  developed  a  shadow¬ 
ing  function  that  can  be  used  with  a  non-Gaussian  probabil¬ 
ity  density  function  [10:788]. 

Approach 

The  approach  of  this  investigation  will  consists  of 
using  Brown’s  shadowing  function 


RIO  ) 


+  {  [q/V  ~  0  P2[q]dq 
V 


I  3 


to  obtain  the  shadowing  function  for  Pyati’s  exponential 
pdf.  In  Eq  (3),  P  is  equal  to  the  cotO,  where  0  is  the 
incidence  angle,  and  P2(q)  is  the  probability  density  of  the 
surface  slope  q.  P^(q)  will  be  derived  for  Pyati’s 


exponential  pdf  using  a  relationship  established  by  Barrick 


between  the  slope  probability  density  function  and  the 
average  of  the  integral  for  the  average  scattered  power. 

The  computed  shadowing  function  will  be  compared  to  the 
results  that  Brown  obtained  for  Barrick’s  'exponential  like’ 
joint  probability  density  function  and  also  with  a  Gaussian 
probability  density  function. 

In  Chapter  II,  an  overview  of  electromagnetic 
scattering  theory  and  an  in-depth  look  at  Brown’s  shadowing 
function  will  be  presented. 

Chapter  III  contains  the  details  of  applying  Brown’s 
shadowing  function  to  Pyati’s  pdf.  Appendices  A  and  B 
contain  the  details  of  obtaining  the  characteristic  function 
of  Pyati’s  exponential  pdf.  Appendix  C  pertains  to  computer 
programs  written  to  numerically  compute  the  shadowing 
f  unc  t ion . 

In  Chapter  IV,  the  results  of  applying  Brown’s 
shadowing  function  to  Pyati’s  exponential  pdf  are  presented. 
Also,  the  development  of  Brown’s  shadowing  function  for  a 
Gaussian  pdf  is  given  in  Appendix  D. 

Chapter  V  contains  conclusions  based  upon  the  theoreti¬ 
cal  analysis  and  computational  results.  Recommendations  for 
future  research  are  also  discussed. 


I I ■  Theor 


This  chapter  contains  the  development  of  Brown's  back- 
scattering  shadowing  function.  To  develop  his  shadowing 
function,  Brown  used  Smith’s  geometrical  shadowing  approach 
and  also,  assumed  that  the  surface  heights  and  slopes  were 
statistically  independent  [10:788].  This  chapter  will 
discuss  the  theory  of  electromagnetic  scattering  first, 
followed  by  the  development  of  Brown’s  shadowing  function. 

Electromagnetic  Scattering 

One  analytical  approach  that  has  been  applied  to  the 
study  of  the  scattering  of  electromagnetic  waves  by  a  random 
rough  surface  is  the  Kirchhoff  approximation  approach. 
Kirchhoff  approximation,  also  known  as  physical  optics, 
approximates  the  fields  at  any  point  on  the  surface  by  the 
fields  that  would  be  present  on  the  tangent  plane  at  that 
point.  This  approximation  is  why  a  large  radius  of  curva¬ 
ture  relative  to  the  incident  wavelength  at  every  point  on 
the  surface  is  required. 

For  the  case  of  backseat  ter ing  ,  when  the  physical 
optics  approximation  is  used  in  the  two  dimensional  case  and 
the  surface  is  assumed  to  be  a  perfect  conductor,  the 
electric  field  at  a  point  P  in  the  far  field  at  a  distance  R 
from  the  surface  is  given  by 


ec>  = 


i  k  R 


'1  HR 


See 


iO) 


e\p  [(  2.)  )  (k  X  +  k  0]  <  l 


k  sin{9)  with  k  being  the  free- 


where  k  =  k  cos(^) ,  k  = 

Z  X 

space  wave  number.  t  is  the  random  surface  height  variable 
and  j  is  the  square  root  of  minus  one  [15:38]. 

The  average  Radar  Cross  Section  ( RCS )  of  the  surface 
that  produced  the  electromagnetic  field  given  by  Eq  (4)  is 
given  as 


'30 


<RCS  > 


2  2 
2k  sec 
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where  X  and  X  are  the  joint  and  marginal  characteristic 
1  2 

functions  of  the  surface  pdf  and  v ^  is  equal  to  2k  sin($). 

Jo  is  a  Bessel  function  of  the  first  kind  of  order  zero. 

A  shadowing  function,  R(3),  is  used  to  correct  the  RCS 
of  the  surface  for  errors  caused  by  shadowing  of  the 
electromagnetic  waves  by  parts  of  the  rough  surface.  The 
RCS  of  a  two  dimensional  surface  that  is  corrected  for 
shadowing  is  given  as 

TC 

<  RCS  >  =  2k2  R (9)  sec2  (  9 ) 

o 

A  more  detail  explanation  of  electromagnetic  scattering 
theory  can  be  found  in  Rosemsn  [15].  Roseman  computes  the 
RCS  of  a  surface  described  by  Pyati’s  exponential  pdf  ami 
then  compares  the  RCS  obtain  with  the  RCS  of  a  surface 
described  by  Barrick’s  'exponential  like’  pdf. 
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Brown ’ s  Shadowing  Function 

Smith  derived  a  shadowing  function  defined  as  the 
probability  that  a  point  on  a  randomly  rough  surface  will 
not  be  shadowed  by  other  parts  of  the  surface.  The  geometry 
of  the  problem  is  shown  in  Figure  1.  The  point  F,  located 
on  the  surface,  is  a  distance  of  £  above  the  mean  plane  of 


Figure  1.  Shadowing  Geometry  of  the  Rough 
Surface.  SF  is  the  Incident 
Ray .  [7:23] 
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the  surface.  The  local  slopes  at  the  point  F  are  pQ  in  the 


x,  z  plane  and  qQ  in  the  y,  z  plane.  The  origin  of  the  .v, 


y,  and  z  coordinates  is  located  on  the  mean  plane  of  the 
surface  directly  below  the  point  F.  Also,  the  incident 
electromagnetic  field  is  orientated  in  the  x  -  0  plane. 
Since  only  the  surface  to  the  right  of  F  can  shadow  F, 
Smith’s  shadowing  function  S(F,$)  is  equivalent  to,  "the 
probability  that  no  part  of  the  surface  to  the  right  of  F 
will  intersect  the  ray  FS"  and  is  given  by  [  19:669] 


S(  F,0  )  =  h  (  F  -  qQ  )  exp 


H 


g(t)dt 


(  7  ) 


In  Eq .  (7),  h()  is  the  unit  step  function,  F  =  cotO 

where  0  is  the  incident  angle  measured  from  vertical,  and  qQ 

is  the  slope  of  the  surface  in  the  y  direction  at  F.  Also, 
gMt)Ar  is,  "the  conditional  probability  that  the  surface  in 
At  will  intersect  the  ray  FS  given  that  it  does  not  in  the 
interval  t”  and  the  result  is  [19:6691 
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Brown  and  Smith  evaluated  Eq  (8)  by  assuming  that  th* 
heights  and  slopes  at.  F,  (  v  =  0),  were  uncorrelated  with 
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those  at  y  -  t.  Brown  also  assumed  that  the  decorrelation 
implied  that  the  heights  and  slopes  were  independent  random 
variables  [10:788-789].  This  assumption  is  valid  for 
Gaussian  distributions;  however,  it  is  not  necessarily  valid 
for  all  other  surface  distributions  [13:136].  However,  this 
assumption  allowed  P^  to  be  written  as  [10:789] 


P3(?  ,  q\  F,  r  )  =  P  (<  )  P2(  q) 


Here,  P  ( i )  is  the  height  probability  density  function,  and 

P  (  q)  is  the  y  slope  probability  density  function.  P  ( q)  is 
2  2 

derived  from  p,  q)  ,  the  joint  pdf  of  the  x  and  y  slopes, 
and  is  given  by  [10:789] 


P2(  q)  ~  P22{p,q)  dp 


Using  Eqs  (8)  through  (10),  Brown  developed  a  shadowing 
function,  R{9),  that  is  defined  as  the,  "probability  that  a 
backsca t te r i ng  specular  point  on  the  surface  will  not  be 
shadowed"  [10:788].  This  function  is  given  by 
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-  1 )  PA  q)  dq 


In  Eq  (11),  t1  -  cot& ,  where  9  is  the  incidence  ingle, 


and  P  A  q )  is  the  probability  density  ,f  t)i,»  surface  si  >pe. 


This  equation  will  be  used  in  the  next  chapter  to  compute 
the  shadowing  function  for  Pyati’s  exponential  pdf, 
Barrick’s  'exponential  like'  pdf  and  a  Gaussian  pdf. 
Shadowing  results  will  be  presented  and  compared  for  the 


III.  Anal: 


Objective 


In  this  chapter,  the  shadowing  function  for  Pyati’s 
exponential  bivariate  pdf,  Eq  (2),  will  be  computed  and 
compared  to  Barrick’s  'exponential  like’  pdf  and  a  Gaussian 


pdf . 


In  order  to  determine  the  shadowing  function  for 


Pyati’s  pdf,  the  joint  probability  density  function  of  the  .v 
and  y  slopes,  ^  (p,  q),  must  first  be  computed  for  Pyati’s 
pdf.  Once  P22  is  known,  it  will  be  used  to  compute  the 
shadowing  function. 


Joint 


Probability  Density'  F  unct  ion 


Barrick  established  the  following  relationship,  valid 
in  the  high  frequency  limit,  between  the  slope  probability 
density  function  and  the  average  of  the  integral  for  the 
average  scattered  power  of  a  rough  surface. 
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In  Eq  (12),  ^(x1»Jr1>  an<^  ^(x2,y2)  are  the  surface 
heights  at  points  1  and  2  on  the  surface.  In  Barrick’s 
notation  £x  =  sin^]  -  sin  C«J  cos  c«j.  £  =  -sin  OO 

sin  C0J  and  K z  =  -cos  -  cos  C^s3  >  where  B  and  B  are 
the  incidence  and  scattered  angles  respectively.  The 
angular  brackets  /  \  denotes  an  ensemble  average.  The  term 
k,  is  the  free-space  wave  number  2w/i.  Also,  the  time 
dependence  exp(-j<*>t)  has  been  assumed  and  omitted  [  2:1  728]. 

The  function  P  on  the  right  side  of  Eq  (12)  is  the 
.joint  probability  density  function  for  the  surface  slopes 
which  is  needed  to  compute  Brown’s  shadowing  function.  To 
evaluate  the  left  side  of  Eq  (12),  the  separation  distance 
between  the  points  on  the  surface  (x^y^  and  (x2,y2>  are 
converted  to  polar  coordinates.  If  r  is  the  separation 
distance  and  <P  is  the  direction  angle  from  (Xj  ,y(  )  to 
( x2 , y 2 )  then 


r  - 


Os  -  x2]2  +  Os  -  y2y 
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x.  - 


x2  =  r  cos<P 


y.  -  y,  -  r  sin'P 


(13) 


Substituting  Eq  (13)  in  the  left  side  of  Eq  (12)  and 
omitting  the  limit,  with  an  understanding  that  it  still 
appLies,  results  in 
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(  1 


Integrating  Eq  (14)  over  <fi  yields 
x> 


2k' 


2  2-. 
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where  v  =  k  t  ,  v  =  k  ^  ,  and  v  =  k  ?  .  Also,  J  is 
x  “xy  y  z  2  o 

a  zero  order  Bessel  function. 

Barrick  showed  that  the  surface  slopes  in  the  x  and  y 
direction  can  be  written  as  ~^x/^2  and  -fv/'z  respectively 
[16:720].  Therefore,  Brown’s  p  and  q  slope  terms  are  aqual 
to  _vx/vz  and  -v^/v  •  Combining  this  with  Eq  (12)  and  (15), 
the  joint  probability  density  for  the  surface  slopes  can  be 
obta i ned . 


P22(P,q) 
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It  Ls  noted  that  the  term  in  the  angular  brackets  is 
equal  to  r  (v^-vj.  Here,  f 2  is  the  characteristic 


function  of  the  joint  pdfs  [5:78-79],  Eq  (16)  becomes 


P22(p,q) 


Equation  (17)  is  the  joint  probability  density  function 
of  the  surface  slopes  that  was  used  to  compute  the  shadowing 
function  for  Pyati's  exponential  pdf.  The  next  section  of 
this  chapter  expands  on  the  characteristic  function. 

Characteristic  Function 

The  second  order  characteristic  function  of  the  random 
variables  <  and  is  defined  as 

x  x 

WV  =  1  I  fVW  exp£  ivi<i+iv2c23  dS  118 

-00  -CD 

where  P_(C  .£_)  is  the  joint  pdf  of  C  and  K  [5:189]. 
z  l  z  i  * 

Also,  <  =  tlx  ,y  ),  K  =  <(x  ,y  )  and  v  =  v  are  used. 

1.  11  Z  i.  Z  Z 

The  joint  pdfs  are  a  functions  of  the  surface-height 
correlation  coefficient  P(r).  Therefore,  the  characteristic 
functions  are  also  a  function  of  the  correlation  coeffi¬ 
cient.  The  variable  r  is  the  separation  distance  between 
the  two  points,  ^  and  on  the  randomly  rough  surface. 

As  r  approaches  infinity,  the  i  n t e rde pendence  of  f  and  ' 
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r  dr  (17) 


l  6 


decreases  . 


are 


At  infinity,  the  terms  <  and  K2 
independent  [3:129]. 

Pyat i 1 s  Exponent ial  Probabi 1 i tv  Dens i ty  Function 
Pyati's  exponential  pdf  is  [14] 


where  I  is  a  modified  Bessel  function  of  the  first  kind  of 
k 

order  k,  h  is  the  mean-square  roughness  and  also,  sgnK  )  is 
equal  to  plus  one  for  f  greater  than  zero  and  equal  to  minus 
one  for  C  less  than  zero.  The  formidable  task  of  solving 
for  the  characteristic  function  is  shown  in  Appendices  A  and 
B  and  the  final  result  is 
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The  joint  probability  density  function  for  the  surface 
slopes  for  Pyati’s  pdf  was  obtained  by  substituting  Eqs  (20) 
and  (21)  into  Eq  (17).  However,  due  to  the  complexity  of 
this  function,  it  could  not  be  integrated  in  a  closed  form. 
Therefore,  the  function  was  evaluated  numerically  on  a  com¬ 
puter,  details  of  which  may  be  found  in  Appendix  0. 

To  gain  confidence  in  t  he  accuracy  of  the  numerical 
integrations,  the  algorithm  was  applied  to  the  Gaussian  pdf 
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for  which  closed  form  solutions  are  known.  Details  of  the 


closed  form  solution  for  the  Gaussian  pdf  are  given  in 
Appendix  D.  The  difference  between  the  numerical  solution 
and  the  closed  form  solution  was  less  than  one  tenth  of  a 
percent.  Therefore,  the  algorithm  was  consider  sufficient 
to  do  the  integrations. 

Correlation  Coef f ic ient 

The  correlation  coefficient,  P(r),  also  called  the 
autocorrelation,  gives  the  correlation  between  K  and  C  . 

One  form  of  the  correlation  coefficient  is  the  Gaussian 
given  by  P(r)  -  exp[-r2/f2],  where  £  is  the  correlation 
length.  Because  of  its  rapid  convergence  and  its  zero  slope 
at  the  origin,  the  Gaussian  correlation  coefficient  has 
found  wide  use  [1:103]. 

Mean-Square  Total  S lope 

Beckmann,  Smith,  Kong  and  others  used  the  following 
equation  for  the  mean-square  total  slope  of  the  rough 
surface  in  their  work  [5:193],  [11:540],  [18:669] 
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and  for  the  Gaussian  case  w  =  2h  /*  .  However,  Bam  ok 

stated  that  since  r2  =  (x  ( -  x  2  +  Cv  t  ~  -v  2  >  the  mean-square 
total  slope  should  be  taken  as  I  1  :  10  4  ] 
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2  2  2 

and  hence,  w  =  4h  /■£  when  P(r)  is  Gaussian.  This 
indicates  a  factor  of  two  variation  between  Barrick’s  and 
Beckmann’s  term  for  the  mean-square  total  slope  of  the  rough 
surface.  This  difference  will  be  discussed  in  the  next 


chapter . 


IV.  Resul ts 


Introduct ion 

In  this  chapter,  following  the  method  used  by  Brown, 
the  shadowing  function  for  Pyati’s  exponential  pdf  and 
Barrick’s  'exponential  like’  pdf  are  compared.  The  first 
section  of  this  chapter  examines  the  impact  of  the  factor  of 
two  variation  between  Barrick’s  and  Beckmann’s  mean-square 
slope.  This  difference  between  Barrick’s  and  Beckmann’s 
mean-square  slope  term  was  discussed  in  the  last  chapter. 
Next,  the  normalized  marginal  slope  density,  P^iq),  is 
graphed  for  Pyati’s,  Barrick’s  and  a  Gaussian  pdf.  Finally, 
the  shadowing  functions  are  graphed  for  all  three  pdfs  for 
several  values  of  the  mean-square  total  slope  term. 


Factor  o f  Two  Difference 

To  avoid  confusion,  Barrick’s  mean-square  surface  slope 

2  2 

term  will  be  denoted  by  w  and  Beckmann’s  by  m  .  Therefore, 
the  following  equation  can  be  written  where  h2  is  the  mean- 
square  roughness  height  and  t  is  the  surface  correlation 
length  . 


I'sirig  this  notation,  Beckmann's  joint,  density  I'unrt  ion 
of  t  h"  surface  slopes  for  the  Gaussian  surface  L  s  |  I  :  tin  | 


exp 


P22(p,q)  = 


2  Tim 


- 

2m  2 

(  2c 


and  that  of  Barrick’s  is  [16:720-721] 
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for 


®  <  p,  q  < 


The  marginal  density,  P  (  q)  ,  is  obtained  by  integrating 


P22<p,q)  over  all  values  of  p.  The  slope  marginal  pdf  for 


Beckmann’s  and  Barrick’s  Gaussian  pdf  are  given  in  Eqs  (27) 
and  (28)  respectively. 
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For  completeness,  the  joint  surface  slope  and  the  sur¬ 
face  slope  marginal  for  Barrick’s  'exponential  like’  pdf  ar- 
presented  in  Eqs  (29)  and  (30)  [10:789). 
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A'  (  )  is  a  modified  Bessel  function  of  the  second  kind  of 
order  one.  The  P2<q)  function  for  Pyati’s  pdf  was  computed 
numerically  from  equations  (17)  and  (20). 

To  demonstrate  the  effect  the  factor  of  two  has  on  the 
results,  both  Beckmann’s  and  Barrick’s  Gaussian  marginal 
slope  densities  are  plotted  in  Figure  2.  The  normalized 
densities  of  Eqs  (27)  and  (28)  are  plotted  as  a  function  of 
P,(q)  times  w  against  the  normalized  slope  q/ w  for  Barrick’s 
and  as  a  function  of  P^  (  q )  ti mes  m  against  the  normalized 
slope  q/m  for  Beckmann’s.  Brown  used  this  method  to  compare 
the  slope  probabilities  of  Eqs  (27)  and  (30)  [3:789). 

It  is  readily  apparent  from  Figure  2  that  the  factor  of 
two  makes  a  difference  in  the  slope  probabilities.  Compar¬ 
ing  Barrick’s  Gaussian  pdf  to  Beckmann's,  Barrick’s  pdf  has 
a  greater  probability  that  a  smaller  q  will  occur  and  a 
decreased  probability  for  a  larger  q.  Small  q  denot.es  a 
sma 1  1  s 1  ope  va 1 ue . 

The  difference  between  the  two  pdf  in  Figure  2  is  gr<  -  it. 
enough  to  re-.»yam  i  ne  Brown’s  results.  Therefore,  the  ne\t 
sec  t  i  .ns  ■  f  this  chapter  will  compare  Barriek's  Gaussian  and 
'exponent  i a  1  Like’  pdfs.  Also,  Pyat  i’s  exponential  pdf  will 
Be  in  -  iaied  in  the  com  par  is  in. 


Slope  P robab 1  I i t y  Dens i ty  Funct  ion 

For  backsea t t e r i ng  to  occur,  the  surface  must  be 
perpendicular  to  the  incident  signal.  Therefore,  the  s  Lope 
q  must  equal  minus  the  tangent  of  the  incident  angle  theta. 
But,  since  /^(q)  is  an  even  function,  q  can  be  written  as 
the  tangent  of  the  incident  angle  when  computing  the 
probability  of  the  surface  slope.  Consequently,  when  theta 
is  equal  to  zero,  the  surface  slope  is  horizontal. 

The  normalized  densities  of  Eqs  (28),  (30)  and  the 

numerical  solution  of  Pyati’s  pdf  are  plotted  in  Figure  3  a 
a  function  of  P(q)  times  w  against  the  normalized  slope 
q/w .  Figure  3  shows  that  all  three  pdfs  have  their  largest 
probability  values  for  small  values  of  q/w  and  that  their 
probability  values  approach  zero  as  q/w  increases.  Small 
values  of  q/w  corresponds  to  small  slope  angles.  Pyati’s 
and  the  Gaussian  pdf  have  about  the  same  probability  for 
small  values  of  q/w,  while  Barrick's  'exponential  like’  has 
a  much  greater  probability. 

Brown  states  that  there  is  an  intuitive  reason  why 
Barrick’s  'exponential  like’  pdf  has  a  much  greater 
probability  for  small  slopes.  According  to  Brown,  having 
the  higher  probability  for  small  slopes  is,  ’  in  agreement 
with  one  intuitive  approach  for  generat  ing  a  surface 
e h a  r a c  t  e  r  i  z e d  by  F. q  (29),  e  .  g  .  ,  one  strongly  filters  all 
surface  height  excursions  below  a  certain  level  to  el  i m  t  na  *  • 

the  po s s i h l  1 i ‘  v  of  large  negative  height  excursions. 


Th  i 


process  increases  the  probability  of  small  slopes  at  the 
expense  of  the  larger  slopes"  [10:789]. 


The  probability  for  large  slopes  should  be  greater  for 
an  exponential  slope  pdf  than  for  a  Gaussian  slope  pdf. 
Looking  at  Figure  3,  Pyati's  exponential  pdf  has  a  slight 
but  significantly  higher  probability  of  larger  slopes  than 
does  Barrick’s  Gaussian  and  'exponential  like’  pdfs. 

However,  the  cumulative  probability  distribution,  area  under 
each  graph,  is  the  same  and  equal  to  a  half  for  each  pdf. 

The  area  under  the  curves  would  be  one  if  they  were  graphed 
from  minus  infinity  to  plus  infinity. 

Figure  4  is  similar  to  Figure  3  except  that  the  mean- 
square  slope  w  has  been  set  equal  to  one  and  the  q  slope 
term  converted  to  its  slope  angle  using  the  fact  that  the 
tangent  of  theta  equals  q.  Therefore,  Figure  4  is  a  graph 
of  the  probability  of  a  slope  angle  verses  the  slope  angle, 
given  that  w  is  equal  to  one.  The  slight  but  s i gn l f i can t 1 y 
higher  probability  of  larger  slope  angles  for  Pyati’s 
exponential  pdf  can  be  seen  in  Figure  4.  The  probability  of 
a  slope  of  70  degrees  with  w  =  I  was  computed  to  be  .0155, 
f >r  Pyati’s  exponential,  .0036  for  Barrick’s  'exponential 
like’  ,  and  .0001  for  the  Gaussian  pdf.  The  next  sent  ion 
:  I  lust  rat  es  Brown’s  shadowing  results  fir  the  three  pdfs 
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Shadow i ng  Function 

This  section  contains  the  results  of  applying  Brown’s 
shadowing  function 


R(0  ) 


0 


P2(  q) 


(31) 


to  various  pdfs.  In  Eq  (31),  /.i  -  cotO  and  (  q )  is  the 

probability  density  of  the  surface  slope.  Brown  defined 
R(£)  as  the  "probability  that  a  backsca t te r i ng  specular 
point  on  the  surface  will  not  be  shadowed"  [10:788]. 

First,  the  results  of  the  shadowing  function  will  be 
compared  when  equations  (27)  and  (28)  are  used  for  (  q )  in 
Eq  (31).  This  will  establish  what  effect  the  factor  of  two 
has  on  the  shadowing  function.  Then,  the  results  of  Brown’s 
shadowing  function  are  compared  for  the  three  different 
pdfs. 

When  Eq  (28)  is  used  for  P  {  q) ,  Eq  (31)  becomes 
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Further  more,  the  backseat  ter ing  shadowing  function  fir 


Beckmann ’ s  Gaussian  pdf,  Eq  < 2~ ) ,  is  given  by 


Appendix  D  gives  the  details  of  the  derivation  of  Eqs  (32) 
and  (33). 


In  Eqs  (32)  and  (33),  h  =cot9,  where  9  is  the 

incident  angle  measured  from  the  vertical  axis.  The 

variables  m  and  w  are  related  by  2m  =  w  =  4h  /<•  ,  where 

2 

i  is  the  surface  correlation  length  and  h  is  the  mean- 
square  roughness  height.  Also,  erfc ( ' )  is  the  complementary 
error  function. 

2  2  2  2 

The  equations  2m  =  w  =  4h  /£  ,  can  be  reduced  to 
h  (m)  =  2h  (w).  Therefore,  the  mean-square  roughness  height 

used  by  Beckmann  is  twice  the  value  used  by  Barrick  when  the 
correlation  length  is  held  constant.  This  greater  height 
implies  a  greater  probability  of  larger  slopes  for 
Beckmann's  pdf.  The  probability  of  the  larger  slopes  was 
seen  in  Figure  2. 

If  the  value  of  the  function  P  ( q)  is  increase  in  an 
interval  in  Eq  (31),  then  the  probability  is  decreased  that 
a  back sea 1 1 e r i ng  specular  point  on  the  surface  will  not  be 
shadowed.  Therefore,  the  shadowing  function  should  be  less 
for  Beckmann's  pdf  because  of  the  higher  probability  of 
larger  angles. 

As  expected,  Figure  5  shows  that  r  he  shadow i ng  funet  ion 
for  Beckmann’s  Gauss  ian  pdf  i s  less  then  Barrick ' s  Gaussian 
pdf.  The  shadowing  functions  for  Beckmann  and  Bar  rick’s 
Gaussian  pdfs,  plotted  in  Figure  ">  with  w  =  m  =  1.11, 


3  0 


.*WV 


^ .* 


Barr  I ck 


28  3a  40  se  60  70  08 

ANGLE  OF  INCIDENCE,  DEGREES 


Figure  5 . 


Backseat  ter  ing  Shadowing  Function  for 
Barrick's  and  Beckmann’s  Gaussian  PDF 
with  W  :  1.41 


shows  that  the  factor  of  two  greatly  influences  Brown’s 


shadowing  function. 

Brown’s  shadowing  function  is  graphed  in  Figure  6  for 
Barrick’s  Gaussian  and  ’exponential  like’  pdfs  and  Pyati’s 
exponential  pdf  with  a  value  of  w  equal  to  .71. 

The  shadowing  function  for  Barrick’s  pdfs,  whose  graphs 
actually  cross  each  other,  have  about  the  same  value.  Where 
as,  Pyati’s  pdf  produces  a  shadowing  function  of  smaller 
magnitude  for  the  same  value  of  w.  The  smaller  shadowing 
function  results  for  Pyati’s  exponential  pdf  are  due  to  the 
pdf's  higher  probability  of  having  larger  slopes  than 
Barrick’s  Gaussian  and  ’exponential  like’  pdfs. 

The  closed  form  solution  for  the  backscatter ing  shadow¬ 
ing  function  of  Barrick’s  'exponential  like’  pdf  was  obtain 
by  Brown  as  [10:790] 

r  x  1  x  f  1  1  1 

R  ( ©  )  =  <  -  K,  (x)  +  -  +  -  K  ( x )  L  ( x )  +  L  (x)  K  (x)  }-  (34) 

[  w  2  2  L  J  J 


The  K  term  is  the  modified  Bessel  of  the  second  kind  of 
n 

order  n.  The  L  symbolizes  the  modified  Struve  functions 

n 

and  the  variable  x  is  given  by  x  =  f/iT/w]  cot.f .  Fq  (01  ) 
was  used  to  graph  Barrick’s  'exponential  like’  pdf  in 


Figure  6 . 


Figure  i  is  a  graph  of  the  shadowing  function  for 

Pyati’s  exponential  probability  density  function  for  severa 

values  of  the  mean-square  slope.  Figure  7  shows  that  the 

2 

larger  the  value  of  the  mean-square  slope  w  ,  the  lower  the 
value  of  the  shadowing  function.  Figure  7  is  used  in 
conjunction  with  equation  6  to  obtain  the  shadow  corrected 
average  radar  cross  section  of  the  surface. 
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Figure  7 


Shadowing  Function  for  Pyati's  Exponenti 
Probability  Density  Function  for  Several 
Values  of  the  Mean-Square  Slope 


V  .  Cone  1 us i on  and  Recommendations 

Conclusions 

The  objective  of  this  study  has  been  to  determine  the 
shadowing  function  for  the  electromagnetic  backscatter  ing 
from  a  randomly  rough  surface  characterized  by  Pyati’s 
exponential  pdf. 

Based  on  a  procedure  developed  by  Brown,  a  shadowing 
function  for  Pyati’s  pdf  was  computed.  The  shadowing 
function  was  found  to  have  a  lower  value  than  the  shadowing 
function  for  Barrick’s  'exponential  like’  and  Gaussian  pdf. 
It  was  determined  that  the  shadowing  function’s  lower  value 
was  due  to  the  high  probability  of  larger  slopes  for  Pyati’s 
pdf.  It  was  also  determined  that  the  shadowing  function  for 
Barrick’s  'exponential  like’  and  Gaussian  pdf  have  about  the 
same  va  lues. 

Recommendat  j  ons 

To  further  the  discipline  of  electromagnetic  scattering 
and  shadowing,  the  following  recommendations  are  made  for 
future  study: 

1.  The  factor  of  two  difference,  between  Barrick’s  and 
Beckmann’s  mean-square  slope  term,  needs  to  be 
reviewed  to  determine  which  author  is  correct . 

This  thesis  showed  that  there  is  a  discernible 
difference  in  the  shadowing  results  when  the 
different  mean -square  f^rms  are  used. 


2.  Further  electromagnetic  scattering  and  shadowing 
research  could  be  pursued  by  building  physical 
surface  models  based  upon  the  three  pdfs  discussed 
and  measuring  their  radar  cross  section  using  a 
laser  set  up  at  different  angles  for  both  the 
momostatic  and  bistatic  case. 

3.  In  this  thesis  much  of  the  computation  was  done 
numerically  on  a  computer  because  of  the  inability 
to  derive  a  closed  form  solution  of  the  character¬ 
istic  function  for  Pyati’s  pdf.  Further  attempts 
should  be  made  for  a  closed  form  solution. 


Appendix  A 


Characteristic  F unc t ion  for  the  Bivariate 
Exponent ial  PDF 

The  purpose  of  this  appendix  is  to  show  how  the 

£ 

characteristic  function  t  {  r  ;p)  for  Pyati’s  two- 
dimensional  exponential  pdf  Eq  (2)  was  derived.  The 
characteristic  function  *  ( v  ,  v  )  of  a  two-dimensional  pdf 
P2(<i,<2>  is  defined  as  [5:189] 

X  x 

VW  =  J  J  P2K!><2)  eXP^  +  JV2r'2^  d<  2  (A_1) 

-  x  -x 

To  find  the  characteristic  function  of  Pyati’s  two- 
dimensional  exponential  pdf,  replace  p2^i>^2>  in  Eq  ( A- 1  ) 
with  Eq  (2)  and  solve.  The  integrals  that  need  to  be 
solved  to  obtain  the  characteristic  function  are  shown  in 


E  q  (  A  -  2  )  . 


In  Eq  ( A  —  2 )  ,  <  and  <2  symbolize  two  random  surface 

heights  which  are  functions  of  .v  and  y.  t  and  r  are 

separated  by  a  distance  r  and  where  P  =  p{  r)  is  the 

surface-height  correlation  function.  I  is  a  modified 

k 

Bessel  function  of  the  first  kind  of  order  k.  Also,  sgn(f  ) 


is  the  signum  function  of  '  . 


f2  <vt,v2;p)  = 


exp 
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At  first  glance,  it  appears  impossible  to  obtain  a 


closed  form  solution  for  Eq  (A-2).  However,  after  much 
trial  and  error,  a  valid  solution  was  obtained. 

The  first  step  in  solving  Eq  (A-2)  is  to  remove  the 
dependence  of  the  integral  upon  the  absolute  value  functions 
and  the  signum  functions.  This  can  be  done  by  using  the 
following  relationships. 


1.2 


-  -  1  for  - v  <  r  <  0 

1.2  1.2 


r  for  0  a  '  ,  ■'  x 

1.2  1.2 


(  A- 3  ) 


s g n ( '  )  =  -  1  for  -a  t  '  ^  <  0 


=  +  1  for  i)  <  ’•  x  2  '  X 


A- 4  ) 


Using  equations  (A-3)  and  (A-4)  the  above  Fourier 
transform  integral,  Eq  (A-2),  can  be  written  as  the  sum  of 
four  integrals. 


t  (  V  ,  V  l  :  I  +  1  a  I  +  1 

’12  1  .1  1  + 


(  A  -  5  ) 


The  values  of  I  through  1  are  given  on  the  next  two 

l  \ 


pages.  A ; so ,  a  substitution  of 
b  r  *  ■  v  i  t  v  . 


h  (  1  -  A’  i  IS  user(  for 


I 


These  integrals  need  to  be  arranged  into  a  form  that 
can  be  integrated  or  for  which  an  integral  is  known.  One 
form  that  works  was  found  in  110:490]  and  is  given  by 
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(  A-  10 ) 


(  A- 1 1 ) 


i  A  -  1  2  ) 


Several  manipulations  must  be  performed  to  get 
integrals  one  through  four  into  the  above  form.  The  first 


step 

is  to  make  a  L  1 

the 

limits  of 

integrat ion 

to 

be 

f  rom 

zero 

to  i n  f l n i t y . 
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four  pages. 


Step  1:  Change  order  of  limits. 
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-  T> 
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Step  2:  Make  all  limits  positive. 
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( A-  1  7  ) 


(  A  -  1  H  ) 


These  integrals  can  be  simplified  by  noting  that  they 


are  really  the  sum  of  two  integrands.  Thus,  if  the 
simplifications  listed  below  are  made,  the  integrals  can  be 
written  in  a  workable  form. 
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A- 2  2  ) 


A- 2 3  ) 


With  the  use  of  these  substitutions  and  equations 
i A  —  1 7 )  through  (A-20),  the  integral  Eq  ( A  —  5  )  becomes 


/  /  D  (A  +  B)  exp  J  C-v/ \  ~vV  drid'' 
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!  A- 2  I  ) 


Grouping  like  terms,  Eq  ( A  —  2 4 }  can  be  rewritten  as 


Tj  =  /  /  D  (A  +  B)  \  exp  j  +v  < 
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( A- 2 5  ) 


Using  Euler’s  identity  exp(j.v)  +  exp(-j.v)  =  2cos(.v), 
Eq  ( A  —  2 5 )  can  be  written  in  the  form 


2  /  /  D  (A  +  B)  cos  (V1/’1  +  V2^2-) 


2  /  /  D  (A  -  B)  cos  [v  <  -  v2C 


(  A  -  2  o  ) 


Equation  (A-26)  with  the  use  of  a  the  trigonometric 


identity  for  the  addition  and  subtraction  of  angles,  which 
states  cosia+l)  =  coset  cos?  x  sinct  sin  3,  can  be  written 
as 


l2  =  2  /  /  D  (A  +  B)  j  cos  (v^)  cos  (v^) 
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{  A  -  2  7  ) 


Expansion  of  Eq  (A-2T)  simplifies  to 
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Now  the  following  substitutions  t  -  C  ( /c  and  n  =  /  c 
are  made  and  the  Bessel  function  is  replaced  with  its  series 
equivalent  Eq  ( A  —  2 9  )  found  in  [20:1  381. 
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Since  v  and  k  are  integers 


T ( i +k+ 1 ) 
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[20:101] 


and  Eq  ( A  —  2 9 )  reduces  to  Eq  ( A  —  3 0  )  . 
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Eq  ( A  —  3 0 )  and  the  A  and  b  substitution  are  used  to 
transform  the  variables  A,  B,  and  D,  given  in  Eq  ( A  —  2 1 ) 
through  ( A  —  2 3 )  ,  into  the  following  equations. 
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Substituting  Eq  ( A  —  3 1  )  through  <  A  —  3 3 )  into  Eq  ( A  —  2 8  > 


Eq  (A-34)  now  has  its  integrals  in  the  form  of 

equations  (A-10)  and  ( A  —  1 1 )  for  which  solutions  are  known. 

£ 

Therefore,  the  characteristic  function  t  for  Pyati’s  two- 
dimensional  exponential  pdf  has  the  following  solution. 
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Where  9  and  0^  are  given  as 
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Finally,  letting  v  =  -v2  =  v  the  requi 
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the  characteris tic  function  i  is  obtained. 
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Appendix  B 


A  C losed  Form  Sol ut ion  for  the  First  Ser les  of  the 
Characteristic  Function 

The  purpose  of  this  appendix  is  to  show  how  a  closed 
form  solution  for  the  first  series  of  Eq  <  A  —  3  5 )  is  derived 
The  first  series  of  Eq  (A-35)  is  repeated  in  this  appendix 
as  Eq  ( B- 1 ) .  The  derivation  of  how  a  closed  form  solution 
for  this  equation  follows. 
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To  obtain  the  variable  P  ,  just  replace  v  with  v,  in 
Eq  (B-2).  The  rest  of  the  equation  is  the  same. 


Using  Euler’s  identity  for  a  cosine,  the  cosine  term, 


cos  ((  k  +  1  >0  )  ,  is  written  as  Eq  ( B  —  3  >  .  The  same  is  done  for 
the  cos((k  +  l)@2]  term  by  replacing  ©  with  0  . 
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(  B-3  ) 


Also,  to  reduce  Eq  (B-l)  to  a  workable  form,  the 
following  substitutions  are  used. 
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The  use  of  equations  (B-3)  through  ( B  —  5  )  allows 
equation  (B-l)  to  be  written  as 


2  , - ,  2  k  r 

1  V_l_  U 


44? 


k  k 


expf.nk  +  1  )  ^  !  *  “\p!  -.j  <  k  1  ;  : 


k  »  o 


j  ex  p  f  j  i  k  ♦  1  !  "  J  >  -  x  p  ! 


-W-l. 


Expanding  Eq  ( B-b  )  by  multiplying  terms  to  remove  the 
brackets  yields 
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Using  the  geometric  series  identity,  )  X  =  ■ _y  which 
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is  valid  when  | X  j  <  l ,  [20:10?!  equation  (B-T)  can  be 

written  as  Eq  ( B-8 )  .  Equation  { B  —  8  >  is  now  in  a  close  form 
since  the  summation  sign  has  been  removed, 
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Euler’s  identities  exp(jx)  =  cos(x)  +  j  sin(x)  and 
exp(-.jx)  =  cos(x)  -  j  sin(x)  are  used  with  the  sub¬ 
stitutions  of  the  terms  !M  =  cos  (6+@  ),  o  =  s  i  n  (  ^  +  9  )  , 

?  =  cos  (  and  $  =  sinl^-0  )  to  obtain  Eq  (  B  —  9  )  . 
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Manipulating  Eq  (B-9)  by  multiplying  each  term’s 
numerator  and  denominator  by  the  conjugate  of  its 
denominator  yields 


(1-P  ) 


2  •'  -  y  -  2p 2 M 2  - 2p  2?<  2 


and 


v. 


o  >  ±7  ±7 

P  -  2--P  P  2lM  +  p  M  +p 


2iP?  -  2p2?  2 -2p2?  2  1 


.2^2  2  4.2  4.2  ! 

a  .h  -  2«fl  +  p  r'  tp  i  J 


Knowing  t n  a  t 

.2  ,2 


c o s  " (  9 l  +  9 } !  +  s  in"i 


os~(9  -9  )  +  sin  (9  -9  )  -  1  allows 
12  12 


( B- 10  ) 


=  1 


F q  ( B- 1 0 )  r  o  be  reduced  to  Eq  (  B- 1  1  )  . 


In  Eq  ( B  —  1 2 )  and  ( B  —  1  3  >  ,  the  trigonometric  identity  for 
the  cosine  of  the  sum  of  two  angles  is  used  to  expand  P1  and 
P  in  terms  of  sine  and  cosine  of  ©.  Also,  Eq  (B-2),  Eq  (B-4) 
and  Eq  { B  —  5 )  are  used  to  evaluate  the  sine  and  cosine  of  ©. 

=  cos  (©  +©  )  -  cos©  cos#  -  sin  ©  sin©  (B-12) 

12  12  12 

P  =  cosi©  -©  )  -  cos 0  cos©  +  sin©  sin©  (  B  —  1 3  ) 

12  12  12 


COS©t 

=  coslcos'l  l/'d  ! 

i  =  1/-" 

( B-14  ) 

COS©2 

=  cosIcosS  1/?  )  ! 

1  =  1/P 

( B- 1 5  ) 

sin©  =  sin  sin1 

l 

sin©?  =  sin  sin1 

Equations  ( B- 1  2  )  through  ( B  —  1  7  )  are  used  to  write  new 
equation  for  P1  and  ?  .  These  terms  are  now  void  of  sire  and 


(  l-P  ) h  v 2 


( 1 -P  ) hv 


( B-  l  T  ) 


cosine  terms. 


1  -  (1-P  )hv  (1-P  )hv2  1  -  (1-P  )  h  v  ( v  2 


1  +  (l-p2)hv  (  1-P 2  ) h v  2  1  +  (  1 -p2  ) 2 h 2 v i v 2 


! B- 1 8 ) 


(B-19  ) 


Substituting  equations  (B-18)  and  (B-19)  into  Eq  ( B  —  1 1 ) 

2  2 

yields  Eq  ( B  —  2 0 )  .  The  values  for  A  B  are  given  in 
Eq  (  B-2 1  )  . 


(1-P  ) 


1  -  (  1-P2 ) 2 h 2  v t  v 2 


2P 2  Q - (  ! -P 2 ) 2 h  2  v  t  v 2 J  +  p4 


2  2  2  2 
1  +  ( 1-P  )  h  vt v2  -  P 

-  2p2Of  l-p2)2h2vtvj  +  P4 


^2B2  =  [l  +  ((l-p“)hvi>2}[l  +  (  (  1 -p“ ) hv2 ) 2] 


( B-20  ) 


(  B-2  1  ) 


Combining  Eq  (B-20)  and  ( B  —  2 l  )  and  factoring  like 
terms  reduces  the  final  equation  to  a  close  form  solution 
for  the  first  series  of  Eq  (  B -  1  . 

if  1  -  !  1  -P  2  )  h  2  v  t  v  ^ 

■l  _  I  ~  2  2  ~  t  2  2  —?.  2  ~~2  2  '  2  ,  2  ~ 

2  [  (1-p  )  h  v  (  v  i  +  h  v  ^  +  h  v  +  2  P  h  v  (  v  ^  +  l 


,  .  2  ,  2 

1  +  (  1  -  P  )  h  v  v 

1  2 

,  2  2  ,  4  2  2  ,  2  2  ,2  2  2  2 

1-P  )  h  v  v.  +  h  v.  +  h  v.  -  Lp  h  V  v.  +  1 


i B- 2  2  ) 


.w 


When  v  =  -  v2  =  v  is  used  in  Eq  ( B  —  2 2 )  ,  the  equation 
can  be  written  as 


2  2  2 

1  +  ( l-p  ) h  v 


A'  = 


,  ,  2  2.4  4  _  ,  2  2  „  2 ,  2  2 

( 1-P  )  h  v  +  2h  v  -  2P  h  v  +  1 


2  2  2 

1  -  (1-P  ) h  V 


2  2,4  4  2  2  „  2  2  2 

( 1 -P  )  h  v  +  2h  v  +  2 P  h  v  +  1 


( B-23  ) 


Equation  (B-22)  and  (B-23)  agrees  with  a  solution 
obtained  independently  by  Roseman  [15).  Also,  when  p  =  0 
Eq  (B-23)  is  equal  to 


A 

a 


l 

V 2  2  r'\"7 

[  h  v  +  1  J 


(  B  -  2  4  ) 


which  is  equal  to  the  marginal  of  Pyati’s  exponential  pdf 
squared  as  it  should  be. 


Compute  X  as  a  Function  of  r 
Store  5000  Data  Points 


The  following  program- was  written  to  compute  the 


shadowing  for  Pyati's  exponential  pdf. 


Variables 

co  r  1  en 
w 
k 
Pi- 
r 
pr 
char 
P 

q 

r  o 


correlation  length 

square  root  of  the  mean-square  slope 
wavenumbe  r 
3 . 141592654 
separation  distance 
correlation  coefficient 
characteristic  function's  value 
surface  slope  x  direction 
surface  slope  y  direction 
shadowing  function 


DEFDBL.  A-Z  '  Double  Precision  used  for  all  Variables 
open  ” c : pam 1 0 . dat "  for  output  as  si 
c  o  r  1  e  n  =  l 


W  =  1.41 

k  r  1  0 

pi  =  4  *  a  t  n  (  1  ) 
v  =  k 

h  =  corlen*w/2 
h  v  =  h  *  v 

dim  qva  1  ( 85  )  , pva  1  ( 8 5 )  , va 1 ue ( 500 5  )  , w ( 5 00 5  ) 
dim  p 2 2(85, 85)  , t e m p ( 8 5 ) 
for  cnt  =  0  to  5000 
r  =  cnt  / 1000 

pr  =  exp(-r*r/(corlen*corlen)  I 
call  chf2(h,v,pr,char) 
w  (  c  n  t  )  -  c  h  a  r 

print  si,  cnt, r, char 
n e x  t  on t 


lose 


/ / / / / /  ll / /  /  /  / / /  /  / 
P,_,ip,q)  and  fy  q) 

/  /  /  /  /  /  /  /  /  /  /  /  ,  /  /  /  /  / 


open  'c  :  pa.)  .  da  t  "  for  output  as  si 
for  qq  =  0  to  80 
q  =  qq/20 

pp  =  qq  1. 1>  8  0 


/  , 


/  20 

for  cnt  =  0  to  5000 
r  =  cnt/1000 
pq  =  sqr ( p*p+q*q ) 
call  bes 1 0 ( k * r * pq , bes ) 
value (cnt)  =  r  *  bes  *  w(cnt) 
next  cnt 

nump  =  5001  :  delta  =  .001 

call  simpsonl value! ) ,nump, delta, result 
p2 2  ( qq , pp )  =  v * v * resu 1 1 / ( 2  * p  i  ) 

p2  2 ( pp , qq )  =  p22(qq,pp) 
print  q , p , p22 ( qq , pp ) 
next  pp 
next  qq 

for  qq  =  0  to  80 

for  pp  :  0  to  80 

pval(pp)  =  p2  2 ( qq , pp ) 
next  pp 

nump  =  81  :delta  =  .05 

call  simpsonf  pval (  )  , nump, delta, result  ) 
qval(qq)  =  2* result 
print  *1  , qq , q , qval ( qq ) 

print  q , qval ( qq ) 
next  q  q 


close 


;/,//, 

1  7 


7  /  /  /  /  /  /////// 7  /  /  /  / 


ttWlTI  T-I^CTK-Ti  R(«  , 


/ / //////// / 


"c;pab.dat"  for  output  as  *1 
for  j  =  1  to  80 

J  J  =  J  *  •  0  5 
for  g  =  j  to  80 
g  g  -  g  *  .  0  5 

cnt  =  1 +  g - j 

t  eir.p  i  on  t  )  =  (  (  gg  /  j  ,i  )  -  1  1  *  q  va  1  i  g  1 

n  e  x  t  g 

npoint  =  80  -  j 
delta  =  .05 

■'a  I  1  simpsonf  f^mpi  )  ,npninf  ,  >1**  I  fa,  resu 
r o  =  1/(1+  result  > 

a  n  g  =  d  (j  -  (  5  "  .  L!  0  5  8  «  a  t  n  :  j  )  )  ) 
p r 1 n '  -  1  ,  ang.ro 

D  i'  1  n  t  .  :  1  ,  ing.ro 


■  v  ■  -  ■  ^  ^  w  »*.  w  ^  ~  ■  *  “  *  ™  k  ’ 


'////// 


Subroutine  that  computes  X. 


sub  chf 2 ( h , v , pr , char )  static 
pi  =  4  *atn (  1 ) 
h  v  =  h  *  v 
pterm  =  l-pr*pr 
c  =  h* (  1 -pr~  2  ) 
termhvp  =  pterm*hv*hv 
ang  =  atn(c*v) 
aterml  =  1  +  ( c*v)~2 

char  =  0 

termll  =  { 1 -termhvp )/( termhvp ~ 2 +2 *( 1 +  p* p )* h v * hv+ 1 ) 

term!2  =  (  1  +  termhvp )/( termhvp ~ 2  +  2  * te rmhvp+ 1  ) 

firstsum  =  ( term 1 1 + te rm 1 2 ) / 2 
fu  =  pr*pi/(4*aterml'1.5) 
sum2  =  f u*s in ( 1 . 5 *ang ) *s i n ( 1 . 5 *ang ) 
n  =  1 

nloopin : 

fu  =  (fu*pr*pr*(2*n+l)*(2*n+lj/(n*(n+L)*4*aterml) 
tempx  =  fu 

t.empy  =  tempx 

if  (tempx  <  3D-9  )  then  goto  nlocpout 

a  =  1 

.)  =  n 
k  =  n  +  2 
1  o o  pin: 

a  2  =  2  *  a 

tempx  =  t  e  m  px  *  (  a  2  -  1  )  *  l  a  2  -  1  i  *  .;  /  i  (a  2+1  )  *  (  a  2  +  1  )  *  k  ) 

tempy  =  t  <-*  m  p  y  +  tempx 

.»  =  .)-  1 

k  -  k  +  1 

l  :,'i+] 

if  <  .)  =  0  or  a  t  1  H  )  t  h^n  goto  Loopout 
g  '.'til  1  ;')pin 
i  "  .pout  : 

■sum 2  :  sum 2  ♦  t  empy  *  s  :  n  (  n  +  1  .  5  i  *  an g  ;  *  s  i  n  t  1  n  +  i  .  5  ;  *  an g  ) 
n  =  it  *  i 
g  ■  •  r  .  i  nl  .pin 

•  \  i  <  •  i  >  >  >  u *  : 

s  e  <  ■  >  ri  •  i  s  u  m  -  H  *  c  *  s  u  m  2  /  <  h  *  p  i  *  p  l  ) 

■•har  =  firs'  sum  +  S'O'nriiisum 


sub  beslO(x,jO)  static 
) 

’  Returns  the  value  of  the  zeroth  ordered  bessel  function 

’Polynomial  approximations  taken  from  Abromowitz  and 
' Stegun  pages  369-370 


j0  =  0  : 

theta  =  0 

:  fO  =  0 

if  (x  < 

-3 )  then 

jl  = 

-99999 

goto  done 

end  if 

dim  j a ( 9 

)  ,  .  j  b  (  9  ) 

if  (x  < 

3  or  x  =  3 

1  then 

goto 

m  1 

else 

goto 

m2 

end  if 

j  a  i  0  )  = 

1  .  00 

j  a  (  1  )  - 

-2 . 2499997 

:  j  a  (  2  )  =  1 

.2656208 

j  a  (  3  )  = 

-  .3163866 

:  j  a  (  4  )  = 

.044447900 

j  a  S  5  )  = 

-  .0039444 

:  .j  a  (  6  !  = 

0002  100 

for  n  - 

0  to  6 

JO  = 

.j  0  +  j  a  (  n  ) 

*(x/3)'(2*n) 

next  n 
goto  d  o  n  e 

m3  : 


a  (  0  )  - 

. 79788456 

jail)  = 

-  .  0000007  7 

.)  a  (  2  )  = 

- 

.00552740 

j  a  (  3  )  = 

-  .  00009  5  1  2 

j  a  (  4  !  - 

.00137237 

j  a  (  5  )  = 

-  .  000  7  2805 

j  a  (  6  )  = 

. 000 14476 

j  b  (  0  )  = 

3 

. oooooooo 

j  b  i  1  )  = 

-  .  7  8  5  3  9  8  16 

.  j  b  (  2  )  = 

- 

.04166 3  9 " 

j  b  (  3  )  = 

-  .  00003954 

j  b  '  4  i  = 

.00262573 

j  b  i  5  .  - 

-.00054 125 

.  j  b  i  6  )  = 

- 

.  001)29833 

.  J  b  1  7  i  r 

.  D'lii  1  3  5  56 

f  i  r  n  = 

0 

f  o  7 

r  h  ^  t 

a 

=  + 

. )  r  .u  n  i  *  (  1 , 

x  )  ■  n  -  1  ) 

r  w<  *  n 

t~‘>r  n  = 

0 

t-. .  >  6 

fo  - 

f  o  +  ]  a  (  ri  )  *  i  I  /  x  )  n 

n  >-■  x  *  n 

J  1  ’<  r  x  ' 

- 

!  /  2  i  *  f  0  * .  in 

i  *  h •• f  a  l 

sub  simpson ( f (  1 ), n , h , result  )  static 
’ Purpose  : 

’  Performs  Simpson’s  rule  of  integration  to  a  function  defined 
’  by  a  table  of  evenly  spaced  values. 

’ Parameters  : 

’  f(  !  -  Array  of  values  of  the  function. 

’  n  -  Number  of  points 
’  h  -  Spacing  between  points 

’  Results  -  Results 
> 

'Check  to  see  if  the  number  of  panels  is  even 
npanel  =  n  -  1 
nhalf  =  i nt ( npane  1  / 2  ) 
n  beg i n  =  1 

result  =  0 

if  npanel  -  2*nhalf  =  0  then 
goto  e  v 
end  if 

’  Number  of  panels  io  odd.  Use  3/8  rule  on  first  three 
’  panels  and  1/3  rule  on  the  rest  of  the  panels. 

result  =  3  *  h / 8  *  (  f (  1 )  +  3  *  f  (  2  )  +  3  *  f ( 3 )  +  f ( 4  j  ) 

nbeg l n  =  4 

e  v  : 

’  apply  1/3  rule,  add  in  first,  second,  and  last  ’.allies, 

result  =  resu  1  t  +  h /  3  *  (  f  (  nbeg  i  n  )  4  *  f  f  n beg  i  n  +  1  )  -t  f  (  n  M 

nbeg in  -  nbeg in  +  2 
if  nbeg i n  =  n  then 
g  o  t  o  ext 
end  if 

'  The  pattern  after  nbegin+2  is  repetitive, 
i }  e  t  nend,  the  place  to  stop, 
nend  =  n-2 

for  l  =  nbeg in  to  nend  s  *  ep  2 

rcsu  1  t  :  re  s  u  It  *■  h  /  3  *  i  2  *  f  1  ;  )  ♦  1  *  f  >  •  ♦  1  • 

ne\t  1 


.  -r.  '  -t 


A  t.  i*.  ,  _  a*  aC  i*.  . V 


Appendix  D 


C 1 osed  Form  S h a d o w i n g  Funct  1  on  for  a  Gaussian 
Probability-  bens  i  t.  y  Function 

This  appendix  will  demonstrate  the  development  of  the 
closed  form  solution  of  Brown’s  Shadow  ms  function  for  a 
Gaussian  pdf.  This  closed  form  solution  of  the  Gaussian  pd  f 
was  used  to  verify  the  numerical  integral  ion  that  was  used 
r  u  compute  the  shadow  i  ntj  function  f  or  Pyafi’s  exponeri  t  i  a  1 
pdf.  The  .joint  Gaussian  probability  density  function,  also 
■  -  a i  1  a  bivariate  Gaussian  density,  is  written  as  i  1  a :  1  o 7  ] 


P . ! 


1 

IF 

i  nh  (1-i.t  J 


- 

>  > 

'  / 1  +  '  ] 

' 

1  12  2 

. 

12  PXP 

2  -  2-.  ! 

L 

(_  'ih  !l-o  j  ! 

_ 

(  n-  i 


The  characteristic  function  of  a  two-dimensional 
pro  t>a  b  l  i  i  t  y  lens  i  t  v  fun  <■  t  ion  is  de  f  i  ri-d  as  |  5  :  1  d  ‘J  | 

+  iv  '  I  d '  a'  i  ! i - i 

2  2J  1  .1 

V  ‘  1  i  rvi  '  no  ha  ra  r  :  s  r  i  ■  f  urc  *  i  ■  n  >  f  •  h**  *  w  - 

me  r:  -  ;  (  n  a  i  j  a  u  s  s  i  a  n  pd  r  ,  replace  !  ■  ,  :  in  >•  u  i  ■  -  C  >  w  i  •  r, 

••  1  i  uni  S' .-  i  \  e  .  The  r  e  1 1  ;  t  of  t  he  :  : :  *  ••  r  v  *  ion  is 


P  (  r  ,  •'  )  exp!  ]  v  ' 

1  i  J  'ii 


r 


i 


V  W'u  r*  v  V  V  •«.  ^  ■%  'N  ' ;  "V  .' 


-r  y,  -  ,  *  ■  *  *  »  ^  ^  y^yw 


Now  let  =  -v  =  v.  Then  i  ( v  , -v  )  can  be  written  as 


I  2  (  v  ,  - v  )  =  exp  [  -h'v'  ^l-pjj 


The  joint  probability  density  f  unct  1  .  >n  o  f  r  >,.•  -.irta  •»* 
slopes  is  obtained  when  the  Gaussian  ■  ■  i .  r r - ! a t  i  n 
coefficient,  P  -  e  x  p  [  -  r  “  / f  “  ]  ,  is  subs  ?  I'a* ;  r,  *  •  j 

which  is  then  placed  into  F  i  l~..  A 1  s  >  •  ,  ?  re-  man.'  i  r.  i  i 


d  e  n  s  1  t.  y  t  ,  w  h  1  c  h  equals  exp!  -  h 


far  t  n«-  ‘ten  -  1  an 


I’orreiat.  ion  <  aoe  f  f  :  c  2  ->n  t  ,  approa  -he  ->  iter-  •  :  n  ♦  he  h  1  tfh 

frequence  limit.  The  joint  pc  f  o  f  the  surf  i  e  i 
written  a  s 


■;  evp  !  -n'V  !_; 


cin  e  on  i  v  r  <•  i  use  »  o  r ro  ,  ■  -,n  t  r  1  hu  • 


ne  e\  p<  iien  *  i  » i  term  exp  i  -  r  .  '  \  an  b* 


t  n i  r :  *  "it 


i  r  s  1  r  w  •  ••  rm: 


"'Mi 


a  ;  |  >  r  \  i  ai  ,a  '  ,  n  ,  r 


•- 


Performing  the  two  integrations,  Eq  ( D-  12 )  becomes 


Ri  0 


1  + 


2  /i 


exc 


-A' 


;A 


e  r  f  c 


A' 

w 


-l 


( D- 1 3  ) 


■hich  reduces  to 


R  i  *  )  = 


1  + 


ex  p 


-// 


erf  c 


! D-  1  4  ) 


Equation  i D  —  i  4  )  is  a  closed  form  solution  to  Brown's 

2  2.2 

snadowing  function  for  a  Gaussian  pdf  with  w  =  4h  /' 

The'*  sha-i  iw ;  ng  function  f  or  Beckmann’s  Gaussian  pdf  can  be 

2  2 

■  no:: hv  subs  *■  i  t  u  f  i  ng  2m"  for  w‘  in  Eq  i  D—  14). 
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